Exact results of a hard-core interacting system with a single impurity 
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The quantum-mechanical problem of a many-particle system with a single impurity in one- 
dimension, interacting by a delta-function, is solved. The wave-function for a bosonic system and 
the related secular equation for the spectrum are obtained. The energy per length is calculated in 
the thermodynamic limit. 
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The quantum-mechanical problem of particles in one- 
dimension with delta-function interaction was studied in 
[1-5] . In ||, the strategy now referred to as 'Bethe ansatz' 
H was successfully applied to this problem. In order to 
obtain the secular equation for the spectrum, periodic 
boundary condition is imposed so that the particles were 
thought of as situated on a circle. Subsequently, some 
other kind of boundary conditions were taken into ac- 
■ count, such as the case of particles being enclosed in an 
infinitely deep potential well Q, the case of with 
modification by addition of an infinitely deep and very 
narrow well at one end of the interval Q or, recently, the 
case of particles being enclosed in a well of finite depth 
H . As a mathematical consideration, had general- 
ized the above problem to the extent within finite Coxeter 
group. 

In the series of a line of studies on exact solutions of the 
above mentioned quantum-mechanical problem, as far as 
we are aware, attention has largely been concentrated on 
the replacement of boundary conditions which are used to 
determine the secular equation. In the present letter we 
consider a system of iV-particles and a single impurity. 
We find the Bethe ansatz wave-function of the system 
and the secular equation for the spectrum. As a result, 
the wave-function is a sum of non-diffractive waves over 
the permutation group Sn and defined piecewisely on 
separate regions. These regions are specified by the per- 
mutation group Sn + i instead of Sn. The energy per 
length of the system is calculated in the thermodynamic 
limit. The problem involved in the present letter differs 
from that of the Kondo Hamiltonian which was ever di- 
agonalized in |n],[^] • 

The Hamiltonian of a system of A^-particles and a sin- 
gle impurity, interacting by a delta-function, reads 
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(i) 

where x' stands for the coordinate of the impurity and 
Xi for that of the ith particle. For simplicity, we start 
discussion with the convention fi := m'/m,h — l,m = 
1. The Hamiltonian ([[]) is invariant under the action of 
whole coordinates translation, so the total momentum of 
the system is conserved. 

The continuity of wave-function and their derivatives 
is determined by the Schrodinger equation on the basis 
of the properties of the potential in the Hamiltonian. In 
order to conveniently obtain the boundary condition aris- 
ing from the delta-function terms in (|l]) , we make a scalar 
transformation x' — > xq = JJix' ' . With this transforma- 
tion, we can combine the first and second terms of (|l|) 
so that they become a standard Laplace operator in an 
(N+l) dimensional euclidean space R w+1 with cartesian 
coordinates (xq, Xi,X2, ■ ■ • , xn). Thus, the Schrodinger 
equation takes the following form 
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2Efe2 +M £ S(Xi-Xj)+vJ2^i- — x ) 
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(2) 

where the summation of the 'kinetic energy' starts from 
I = instead of I = 1. 

After removing the delta-function terms from the left 
hand side to the right hand side of (|^), one takes inte- 
gration of (|^) over a 'Gauss box' that is cut into halves 
by either the hyperplane {xi — Xj = 0| i = 1, 2, • • • , N} 
or {xi - xo/^/JI = Q\i = 1, 2, ■ • • , N}. With the help 
of Gauss integral theorem, we obtain the discontinuity 
relations for the derivatives of wave-function along the 
normal of a hyperplane P a = {x\ (a\x) = 0}, where a 
stands for the normal vector and (a\x) stands for a scalar 
product of two vectors a and x = [xq, X\, x%, ■ ■ ■ , xn) i.e. 

lim {a ■ V[^>{xr a \ + ea) - ip(x( a -\ — ea)]} = 2ci/)(X( a )) 

(3) 

where c = u for a — e; — e ? - and c — v for a — et eo 

N d 

(i = 1, 2, • • ■ , N); V := ^2^1 — ; X( a ) E P a and {e t } is 

the standard orthogonal basis. 

It is convenient to discuss the present problem in 

the frame of reference of the impurity. This requires 

a coordinate transformation (xq, xi, X2, ■ • • , xn) h ~ > 

(2/0,2/1,2/2, ••• ,Un) given by y Q := x' = ^ x o, and 
Ui := Xi — x 1 = xi — Xq. It results that 



ip(x ,xi,x 2 , ■ ■ ■ ,X N ) = Elp(x ,Xl,X2, ■ ■ ■ ,x N ) 
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— = —{t. V t^-) and — = — . Obviously, 

dx y/Ji oyo t-J dyt dxi dyi 

d 9 9.., „ 

— — = > — h — — is just the operator of total mo- 

oy j^i dxi Ox' 

mentum (apart from a factor — i). In terms of these coor- 
dinates (yo, yi, 2/2, • • • , Vn), the Hamiltonian (jl|) and the 
discontinuity relation (||) become 

I N 8 2 \ ( d N N N 

» = -i § m - % (js - § % j + »,£ '<» - *> + 

(4) 
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e _3+ I % ~ 9^J ^ V(o) + ^ ~ ^ (?/(Q) ~ = 2U ^ V ^ 

(5) 



and 



lim 

c-0- 
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where a stands for the normal vector of the hyperplane 
{y\Ui ~~ Vj = 0} and a' for that of the hyperplane 
{y\Vi = 0}- Clearly, (||) is invariant under the trans- 
lation of j/o ~ > 2/0 + c or under any permutation of the 
coordinates j/j for i = 1, 2, • • • , JV. Because of the trans- 
lational invariance of j/o, we may set ip(yo, yi, • • • , j/at) = 
e lKyo ip(yi, • • • , 2/at), where K is a constant, i.e. the total 
momentum is conserved. Then the Schrodinger equation 
is reduced to a differential equation for <f(yi, yi, • • • , J/jv)- 
On the domain R N \{Pp} it becomes 



N q2 1 w a 

9^ 1 9 



2^92/f 2 M ^0» 



V(Vi,y2, ■■■,Vn) = Eip(yi,y 2 , J/jv), 
(7) 



where { IP^} stands for the set of hyperplanes {y\ yi—yj = 
0} and {y\yi — 0} for i,j = 1,2, • • • , AT. Obviously 
these hyperplanes partition R w into finitely many re- 
gions and a wave-function must be piecewisely contin- 
uous function defined respectively on separate regions. 
The regions into which the hyperplanes {{y\ yt — yj = 
0}\i, j — 1, 2, • ■ ■ , N}} partition R N can be specified by 
elements of a permutation group Sn- In other words, 
different regions have different orders for the coordinates 
(j/i) 2/2, • • • , 2/ at) if their components are arranged from 
large value to small ones, e.g. j/i > 2/2 > 2/3 > • ■ ■ > 2/jv, 
2/2 > 2/i > 2/3 > • • • > Vn etc. Moreover, each of such 
regions is partitioned into N + 1 separate regions by the 
other type of hyperplanes {{y \ yi = 0}}, For example, the 



3 



region y\ > y2 > ■ ■ ■ > Un, is partitioned into the regions 
> yi > y 2 > ■ ■ ■ > un, Hi > > y 2 > ■ ■ ■ > Un, ■ ■ • and 
Ui > 2/2 > ■ ■ ■ > yN > 0. These regions can be specified 
by a permutation group Sn+i- It is easy to find a spe- 
cial solution for ([?]), i.e. a plane wave solution y>{k}(y) ~ 
e *(*l») i here (k\y) = kiyi + k 2 y 2 H h k N y N . The cor- 

j JV j JV 

responding energy is E 1 = - /c 2 + — (^J fcj — if) 2 or 

i=l " i=l 

more neatly, 

where a parameter X = K — ^ i s introduced. 

Because of the permutational invariance of the Hamil- 
tonian (0), we look for solutions of the Bethe ansatz form: 

<p T (y) = E ^(^e^M (9) 

a£S N 

where ak stands for the image of a given k := 
(hi, k 2 , ■ ■ ■ , kff) by a permutation a S 5jv and the co- 
efficients A(c, t) are functionals on Sn ® Sn+i- It is 
worthwhile to mention that the summation in (^) is taken 
over the permutation group Sn whereas the various re- 
gions on which the wave-function is defined are specified 
by elements of Sn+i- This is different from the Bcthc- 
Yang ansatz [|. As far as the Bethe ansatz solution (§) 
is concerned, the boundary condition (^) becomes 

d A 

( ~ — i-) y(v)\ yi ^o+ - <p(y)\ Vi ^o-] = 2, "ip(y)\vi^o- 

oyi fj- 

(10) 

For a bosonic system, the wave-function is supposed to 
be symmetric under any permutation of the coordinates. 
Because any element of permutation group Sn can be ex- 
pressed as a product of the neighboring interchanges ac 
(• • -iVhVi+i, •••)!->•(•• ■ ,y t+1 ,y t ,- ■ •), i = 1, 2, • • • , N- 1, 
it follows that {<Jiy){y) = <p(y)- Since <p is a scalar func- 
tion, ((Ji(p) is well defined by f(a~ 1 y). Thus both sides 
of the equation can be written out by using (^|). In terms 
of the evident identity (ak\a~ 1 y) = (er,;<7fc|y) and the 
rearrangement theorem of group theory, we obtain the 
following consequence: 

A(a,aiT) = A(aia,T), fori = 1,2,- • ■ , AT - 1. (11) 

Substituting (||) into the boundary condition (|l0]) and 
using the continuity condition for the wave-function on 
the hyperplanes, we obtain the following relation (when 
yi is nearest to zero among all {yj,j = 1, 2, • • • , N}) 

x v — i\(ak)i — \/u] ,, . , „. 

A(a, a N r = . 1 .yU a, r , 12 
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where erjv € Sjv+i but ^ Sn- The relations (^Tj) and 
( |l2| ) relate A-coefficients between different regions, so we 
direct our attention to one of the regions, e.g. y\ > yi > 
■■■y N >0. 

Because the region with • • ■ yi < yi+\ • • • is next to the 
region with • • • y i+1 < y, • • • for yiy i+ \ > (• • • y t < < 
iji+i • • • is not next to the region • • • yi+\ < < •), 
the boundary condition (||) and the continuity condition 
for the wave-function on the related hyperplanes give a 
relation between ip T (y) and fcnriy) (i = 1) 2, • • • , N — 1). 
After writing out the relation in detail we obtain the 
following: 

i[{ak)i - (ak) i+ i][A(a,T) - A(aia,r) ~ A(a,air) + Afaa, (7^)} = 2u[A(a,r) + A(cr j; cr, r)]. 

(13) 

With the help of (O) this relation gives 



A{(na,T) = -Yi(ak)A(a,r), 
= u- l [(ak)-{ak) +l ] 
V ' u + i[(ak)i-(ak) i+1 ] V ' 

As a result, all the ^-coefficients are determined up to an 
overall scalar factor by using (J14] ) repeatedly. The con- 
sistency of successive use of (|14|) is guaranteed because 
the Y's in ( |l4|) satisfy Yang-Baxter equation. 

In the above we have obtained the Bethe ansatz so- 
lution of Schrodinger equation on the basis of bound- 
ary condition arising from delta-functions in Hamilto- 
nian. Now we turn to determine the secular equation 
of the spectrum {k}. Suppose the system is situated on 
a circle of length L, the wave-function must obey the 
periodic boundary condition ip(yi, ■ ■ ■ , — L, ■ ■ ■ , y^) = 

<p(yi,-'-,Vi,-'-,VN)- if y = (yi, - ■ ■ ,yi, - ■ ■ ,vn) is a 

point in the region specified by r S SV+i, as a result 
of periodicity, y' = (yi, ■ ■ ■ , y» — L, ■ ■ ■ , i/jv) must be a 
point on another region specified by 7T € Sn+i, where 
7 = ctjvA and A := <tn-i ■ ■ -02(71 G Sjy Then the re- 
quirement of the periodic boundary condition is explicitly 
written as <p T {y) — ^71- (?/')• After writing this relation 
in terms of the Bethe ansatz solution (^|), we find that 
the periodic boundary condition is guaranteed as long as 
the relation A(a, 'yr)e~ 4 (' Tfe )i£ = A(a, r) is imposed. This 
implies the following equation after using (pT|), ( |l2| ) and 
( |l4| ) repeatedly, 

l (<rfc) 1 L = / r] N v - i K< Tk )i ~ A /M] TT u-i[(ak) 1 - (o-fc)j] 
1 ' v + iKak^-X/^j^u + iKak^-iak^] 

(15) 

As ( p"5| ) holds for any <r G S^r, (<rfc)i will take 
k±, k2, ■ ■ ■ > fcjv respectively. So (15) represents TV distinct 
equations. 

Taking the logarithm of ( |l5| ) , we have a system of cou- 
pled transcendental equations, 
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2tt, 
L ' 



2 
L 



tan 



X/p 



N 

tan" 



(16) 



where the choice of range for tan -1 is (— 7r,7r), and In- 
takes integer values or half integer values accordingly as 
N is even or odd. ( |l6| ) is the secular equation for the 
spectrum and the {Ij} play the role of the quantum num- 
bers. As a result of periodicity, the total momentum is 
quantized i.e. K — n(^r), n is an integer. In particular, 
when ji = 1, the secular equation recover the known re- 
sult H of an N + 1 particles . From (||) we learned that 
eigenstates related to {A;i, fe, • • • , A;jv} and those related 
to {A, &2, • ■ • , fcjv}, A, • • • , fcjv} etc. have the same en- 
ergy if /i = 1. However, the (N + 1) fold degeneracy is 
broken by the replacement of a particle by an impurity. 

Although the transcendental equation is difficult to 
solve, we are able to do some meaningful calculations 
for N — > oo, i.e. thermodynamic limit. Suppose that 
{fci, fca, • • • , fejv} i s a self-consistent solution set of (16). 
Consider the function 



I(k) 



1 

2^ 



kL + 2 



tan 



k — A / ' fi 



N 

tan" 



(17) 



which is a monotonously increasing function of k. 
Clearly, when I(k) passes through one of the quantum 
numbers Ij, the corresponding k is equal to kj, one of 
the roots. However, there may exist some integer (or 
half integer) values for which the corresponding k is not 
in the set {fcj., fca, • • • , fcjv}, then such a k is called a 
'hole'. A smooth, positive-definite density (per length) 
describing the distribution of roots and holes is natu- 

1 dl{k) 

L dk 

we get an expression for p(k). In the thermodynamic 
limit (i.e. sufficiently large N) one can make a replace- 
ment, namely, limjv_ >00 J2i=i f{h) = J dkp(k)f(k) - 
Sj=i f{hj)> where hi, h%, ■ ■ ■ , h m are positions of holes. 
Then we have 



rally defined as p(k) 



Differentiating (ft 



v 2 + (k- \/pf 



E 



(k 



(18) 



This can be solved by means of the technique of Fourier 
transformation. After some calculations, one gets the 
density of roots and holes as a Fourier integral. Then 
the energy per length of the system is obtained as 



E = 



X 2 



1 



2p,L 4tt 



fc 2 e -i£fc 

I - fi -«l£l 



(19) 



dk'p(k')- 



{k-k'f 



d£dk. 
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for repulsive case u,v > 0. 
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